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Recent investigations have highlighted the failure of a sum of terms to represent the dispersion 
interaction in parallel metallic, anisotropic, linear or planar nanostructures [J. F. Dobson, A. White, 
and A. Rubio, Phys. Rev. Lett. 96, 073201 (2006) and references therein]. By applying a simple 
coupled plasmon approach and using electron hydrodynamics, we numerically evaluate the disper- 
sion (non-contact van der Waals) interaction between two conducting wires in a collinear pointing 
configuration. This case is compared to that of two insulating wires in an identical geometry, where 
the dispersion interaction is modelled both within a pairwise summation framework, and by adding 
a pinning potential to our theory leading to a standard oscillator-type model of insulating dielectric 
behavior. Our results provide a further example of enhanced dispersion interaction between two 
conducting nanosystems compared to the case of two insulating ones. Unlike our previous work, 
this calculation explores a region of relatively close coupling where, although the electronic clouds 
do not overlap, we are still far from the asymptotic region where a single power law describes the 
dispersion energy. We find that strong differences in dispersion attraction between metallic and 
semiconducting / insulating cases persist into this non-asymptotic region. While our theory will 
need to be supplemented with additional short-ranged terms when the electronic clouds overlap, 
it does not suffer from the short-distance divergence exhibited by purely asymptotic theories, and 
gives a natural saturation of the dispersion energy as the wires come into contact. 
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I. INTRODUCTION 

The dispersion interaction constitutes the outer, at- 
tractive part of the van der Waals (vdW) force (in 
the non-retarded regime, for charge-neutral, non-polar 
species). It can be explained via the interaction between 
small instantaneous dipoles arising due to a mutual dy- 
namic polarization of the electron clouds (see e.g. [1]). 
These attractive forces, although weaker in magnitude 
than ionic or covalent bonds between atoms or molecules, 
are ubiquitous across nature and can play a central role 
in soft matter and biophysical systems. Examples range 
from protein folding^i^i^ to the adhesive properties of 
gecko feeli^i^. While the present work cannot describe 
the contact regime where the van der Waals force is 
strongest, it does suggest that noncontacting regions may 
contribute more than previously suspected to the ener- 
getics, in cases where pi-conjugation, for example, may 
lead to near-metallic behavior in low-dimensional struc- 
tures. 

Dispersion forces also play an important role in the 
rapidly progressing area of nanoscience. Nanotubes are 
quasi-one dimensional structures, with electronic prop- 
erties determined by their composition. The metallic or 
semi-conducting behavior depends on the hclicity of car- 
bon nanotubes. Nanotubes composed of boron nitride 
are primarily semi-conducting, with a wide band gap and 
stable electronic properties^'*. Armchair (n,n) carbon 
nanotubes are metallic. The vdW attraction facilitates 
the self-assembly of single-walled nanotubes into bundles 



or ropesS. and being the primary inter-tube attraction is 
important in a description of the cohesive energetics of 
nanotube bundles^-. 

Recently the differing asymptotic behavior of vdW 
interactions for conducting systems, compared to insu- 
lating systems with identical geometry, has been dis- 
played for various systemsiiiiii^ii^. The simplest stan- 
dard techniques applied in computing the vdW inter- 
action are based on the pairwise summation of all 
contributions^ between microscopic elements sep- 
arated by distance R, and so inherently assume local 
properties of the two materials. More sophisticated 
techniquesi^*^ when coupled with the usual assump- 
tion of a local bulk-like dielectric function within the 
boundaries of the constituents, still lead to an inter- 
action asymptotically equivalent to a sum of con- 
tributions. Recently, progress has been made in for- 
malisms that describe both the contact region of elec- 
tron cloud overlap, as well as a dispersion contribution at 
large distanceaiiii^ii^i^i2iiS^i2^i2ii2^. However these the- 
ories predict the usual asymptotic Y,R~^ behavior and 
so could still be improved with respect to description of 
low-dimensional systems with a small or zero electronic 
energy gap. 

For a determination of the vdW interaction between 
highly anisotropic systems of metallic nature, a pairwise 
summation method fails to incorporate the long wave- 
length, incompletely screened electron density fluctua- 
tions. These lead to unexpected power laws for the vdW 
interaction as a function of separation--. The effects of 
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this physics appear in calculations on parallel ID electron 
gases^^, on two dimensional electron gases^^'^'*'^^, and in 
graphene sheets^^. More recently this has been high- 
lighted as a more general phenomenon^^. The greatest 
underestimation of the vdW interaction occurs between 
distant parallel quasi one-dimensional conductors— . 

The coupled plasmon approach addresses the short- 
falls of the pairwise summation technique, incorporating 
polarization between multiplets of atoms in the system 
and also electron movement within the wires. In mod- 
elling the electronic motions as coupled plasmons, large 
electron displacement across many atoms in a system is 
allowed to occur, providing a less constrained represen- 
tation of electron movement. Plasmon frequencies can 
be determined by application of electron hydrodynam- 
ics and the vdW dispersion energy then corresponds to 
the separation dependent part of the sum of zero point 
plasmon energies. This approach has been employed pre- 
viously in determination of the vdW interaction between 
systems of thin parallel metallic plates and also thin par- 
allel metallic wires. The sum of plasmon energies is also 
known to provide an approximation to the full correla- 
tion energy in the Random Phase Approximation (RPA) . 
(See e.g. [13]). A full RPA energy calculation would thus 
provide a natural and seamless extension to the present 
method, valid at all separations including the regime of 
full electronic cloud overlap. 

For two thin parallel metallic plates of infinite area, 
separated by a distance greater than the thickness of 
the plates, Bostrom and Sernelius^^, and Dobson and 
Wang2^ have applied the coupled plasmon approach to 
obtain a dispersion energy of dependence D~^l^ on the 
separation distance D of the metal slabs. This result is in 
agrement with the appropriate limit of microscopic ran- 
dom phase approximation calculations for a pair of two 
dimensional electron gases performed by Sernelius and 
Bjork^^. A pairwise summation analysis for this system 
of thin parallel metallic plates gives a dependence 
of the vdW attraction, notably smaller at large separa- 
tion than the D~^l'^ dependence derived by applying the 
coupled plasmon approach. An equally unconventional 
result was obtained for the interaction between graphene 
planesiiii^,. 

Another example system in which the pairwise summa- 
tion analysis was shown to underestimate the dispersion 
interaction for conductors is that of two parallel wires 
of infinite lengt h^^'^^i'^°i^^ . For a pair of parallel, con- 
ducting, infinitely long wires, the dispersion interaction, 
calculated by a standard sum of R^^ contributions be- 
tween microscopic elements separated by distance i?, has 
a r)~^ dependence on the separation distance, D, of the 
two wires. However an analysis of the zero-point en- 
ergy of the delocalized coupled one-dimensional plasmon 
modes parallel to the long axis, modelled for wires of 
length 1/ >> D, finds the vdW interaction to have a 

D-^ {ln{D/A)-^/^^ (1) 

dependence on separation distance, D. A denotes the 



smearing radius of a wire, representing the finite extent of 
the electronic wave-function or electronic density fluctua- 
tion on the wire in the direction perpendicular to the long 
axis. The analysis assumes an electron mean free path 
greater than the separation distance of the wires, which 
can be satisfied by bismuth nanowires'^^ and conducting 
nanotube O"" ' !^'^ . The vdW interaction ^ is almost three 
powers of D greater than that obtained in a standard 
pairwise summation framework. This highlights the im- 
portance of including electron movement and screening 
from subsequent polarization of electron pairs in conduct- 
ing systems. The occurrence of an enhanced dispersion 
interaction in metallic systems might have repercussions, 
particularly in nanotechnology. The above considerations 
suggest metallic and non-metallic nanotubes could expe- 
rience different cohesive forces, for example though anal- 
ysis beyond the present methods would be needed to ex- 
plore this possibility in the case of electron cloud overlap. 
In the present paper we explore another casc'^^ where a 
sum of R^^ contributions does not yield a correct de- 
scription of the vdW interactions. We consider two lin- 
ear structures ("wires"), each of length L, in a coUinear 
"pointing" configuration as shown in figure 1. D is the 
separation between the near ends of the wires, while Dcm 
is the distance between the centers of the wires. A is an 
effective radius of the wire, discussed further below. 

Our motivation for looking at this case was threefold. 
Firstly, there is intrinsic interest in the cohesive prop- 
erties of nanostructures of all kinds, and the present 
work is particularly relevant to the interactions of car- 
bon nanotubes, a technologically important case. Sec- 
ondly, recent work^^ has suggested that van der Waals 
forces may be involved in the surprisingly strong force 
that tends to make the tips of iron micro-whiskers grow 
towards one another during fluidized-bed iron ore reduc- 
tion processes. While it is not clear that the electron 
mean free path in these cases is sufficiently long to vali- 
date the assumptions of the present work, it is interesting 
that our approach can predict vdW forces between linear 
conducting systems in the appropriate "pointing" geom- 
etry, that are enhanced compared with standard vdW 
theory. Thirdly, a numerical investigation of wires of 
finite length allows us to look at a region of relatively 
close coupling where, although the electronic clouds are 
not permitted to overlap, we are still far from the truly 
asymptotic regime where a single power law describes the 
force. We find that the dispersion interaction between the 
wires shows strong differences between metallic and insu- 
lating cases even in this non-asymptotic regime, reinforc- 
ing asymptotic resultsiiii^'^^'^iiS^ already known in other 
geometries. We further find that the greater interaction 
at large distances in the metallic case does not imply a 
lesser interaction at small distances in the metallic case, 
contrary to a false expectation based on a single power 
law at all distances. 

The paper is organized as follows. The analytic work 
will be presented in section [III and section HITl then out- 
lines the numerics. Our results and discussion compare 
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FIG. 1: 



the dispersion interaction between wires in the "point- 
ing" geometry for insulators and conductors and are pre- 
sented in section IIVI The principal findings and future 
prospects are discussed in the final section. 



II. ANALYTICS 

We evaluate the dispersion interaction between the two 
wires by the coupled-plasmon approach (see e.g. Ref. 
[12], [13], [14]). Plasmons are quantized versions of the 
wavelike collective motions of the electrons. The equilib- 
rium /ground state of the wire is described by a constant 
electron number density per unit length, uq. Plasma 
waves entail a perturbation 6n(X,t) to the number den- 
sity, where we have defined a variable X (see Fig. [T|) to 
label positions along the wire. 

Here we neglect electronic radial and azimuthal mo- 
tions - i.e. those perpendicular to the long axis. These are 
frozen out by quantum effects in atomically thin systems, 
but can be present in other cases. We do not consider 
them here because, for the thin linear systems of interest 
here, they consist of spatially constrained electronic mo- 
tions leading to conventional van der Waals force laws. 
By contrast the unconstrained electronic motions along 
the wire will be shown below to lead to an unconven- 
tional, enhanced long ranged van der Waals interaction 
that is the principal focus of attention here. For example, 
the azimuthal electronic motions have been considered by 
Rotkin and Hess^, in the case of two long parallel nan- 
otubes. They considered only small separations between 
the tubes and found these modes to give a more rapidly 
decaying van der Waals interaction as a function of nan- 
otube separation D, compared to the purely longitudinal 
modes investigated in [11]. Thus we expect that our ne- 
glect of the electronic motions perpendicular to the long 
axis will not affect the long-ranged forces that we are 
investigating here. 

To describe the electronic motions along the wire we 
use a conventional hydrodynamical picture, valid even 



for degenerate quantal electrons in the limit of long- 
wavelength perturbations. During the plasmon motion, 
the electron fluid element that was at position X (see 
Fig. [T|) in the unperturbed state of the wire (bold lines 
in Fig. 1) is displaced to position X + R{X,t) (dotted 
lines in Fig. [T]). In the simplest hydrodynamical model 
(correct to lowest order in the wavenumber Q of the 
wave) the motion of the fluid element is described by 
Newton's second law for a a free mass under the action 
of a mean potential energy function ^{X,t) generated 
by the Coulomb interaction with other fluid elements: 
M(PR{X,t)/dt^ = ~d'i>{X,t)/dX = F. Here M is the 
electron effective mass for motion along the wire. 

An element of electron fluid at a point on any wire will 
experience a Coulomb potential due to electrons at every 
other point on that wire and also due to those at every 
point on the adjacent coUinear wire. Considering plas- 
mons in one dimension, evidently the Fourier transform 
of the Coulomb potential does not converge. However, 
by recognizing that the one electron wavefunctions are of 
finite extent in the direction perpendicular to the wire, 
we can examine a symmetrically 'smeared' version of the 
wire with finite smearing radius, A. While we apply a 
radially smeared pair potential, expressed explicitly as 



m-x'\) = 



((X-X')2+A2) 



1/2 



(2) 



in our work to follow, any expression that saturates the 
Coulomb interaction for |r| < ^ would give similar results 
for the long- wavelength fiuctuation phenomena that drive 
the unusual effects to be explored here. 

The force on an element of electron fluid at a point 
in time can be expressed using the radially smeared pair 
potential, 0, as: 



M 



cpR{X,t) _ d 



-D/2 



D/2+L 



Sn{X',t)ij){\X - X'\)dX' 



D/2-L 

Sn{X',t)(t>{\X - X'\)dX' 



D/2 



(3) 



This applies for either D/2 X D/2 + L or -{D/2 + 
L) ^ X ^ —D/2. We apply a linearized form of the 
continuity equation, 

5n{X,t) = -^{n,{X)R{X,t)), (4) 

to describe the perturbation to the electron number den- 
sity in terms of the equilibrium number of electrons per 
unit length of wire, no{X). 

Plasmon movement within metallic wires is conflned 
by requiring zero electron displacement at the wire ends. 
An alternative analysis, in which the zero boundary con- 
dition occurring at the wire ends is naturally inherent. 
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could be obtained by applying a sine-basis decomposi- 
tion. Note that different boundary conditions, such as 
allowing the electron gas to move out over the ends of 
the uniform positive background, which would then cause 
a restoring force, could be applied instead, and possi- 
bly would provide a further alternative description of the 
plasmon modes. 

We seek time-periodic separable solutions of ^ in the 
form 

R{X, t) = R{X) expi-iflt) . (5) 

The left-right symmetry evident in Fig.[l]requires that 
there exist even solutions for which R{X) = R{—X) and 
also odd solutions for which R{X) — —R{—X). We used 
this property and put ([5]) and (|4]) into ([3]). Integrating by 
parts and using the explicit form of the smeared Coulomb 
potential ([2]) we obtained 

M^-'RiX) _ /■^^^^'^.^/^ / -2(X-X'f + A^ 
noe^ ^ Jd/2 \{{X - X'f + A^f^^ 

{{X + x'y + A^y^^ J 

(6) 

Eq. © holds for either -{L + D/2) < X < -D /2 or 
D /2 < X < L + D/2, and is an eigenvalue equation for 
the frequencies fi^ of self-sustaining plasma oscillations, 
(plasmons) corresponding to coupled charge density fluc- 
tuations on the two wires. 



A. Evaluation of the dispersion energy 

The dispersion energy E^^w is then the part of the to- 
tal plasmon zero-point energy that depends on the sepa- 
ration D between the wires. Thus 

i^.dw(i5) = ^^(r!,(i^)-l],(D^oo)) . (7) 

3 

Application of numerical techniques is necessary to 
solve for the eigenfrequencies fi^ , as the broken trans- 
lational symmetry of the coUinear wire geometry with a 
gap prevents Fourier transformation being used to solve 
the integral equation ^ analytically in k-space. 

B. Incorporating Internal Pressure 

Although the eigenvalue problem ^ derived from the 
expression for the force felt on an element of electron 
fluid at a point could now be solved to give us odd and 
even solutions for the plasmon frequency, we extend this 
to incorporate the contribution to the force felt by an el- 
ement of electron fluid from the internal pressure of the 



electron fluid. This electronic pressure is caused mainly 
by velocity deviation of electrons in a fluid element from 
the average velocity. It is the electron pressure term that 
distinguishes our formalism from a purely classical New- 
ton II approach, and it introduces the quantal and Pauli- 
principle physics (electron degeneracy pressure) into the 
formalism. We expect it should give higher frequen- 
cies to the plasmon modes with rapid spatial variation. 
Without the pressure term the theory correctly describes 
the response of a free electron gas in the limit of long 
wavelength, and it agrees with the corresponding limit 
of the quantal Lindhard response. This limit was ad- 
equate in [12] because there the treatment was limited 
to the asymptotic regime of distant interactions between 
infinitely long wires, where the long- wavelength fluctu- 
ations completely dominate. Here we want to explore 
wires of flnite length at finite separations. While the 
new phenomena which we expose will still turn out to 
be due to long-wavelength fluctuations, a detailed treat- 
ment needs to account, for example, for the reduced po- 
larizability of a short wire of length L resulting from the 
requirement to excite wavelengths of 0{L~^) in order to 
polarize its electron gas. For this reason we introduce 
the plasma pressure, in the simplest available model. In 
some regimes it also aids in achieving numerical stability 
of the solutions. 

The excess pressure due to a density perturbation 6p 
is usually modelled from an analysis of the free Fermi 
gas^SiSl^iSQ^^ as SP = MB^Sp. Here Sp is the three 
dimensional density perturbation and _B is a velocity of 
the order of the Fermi velocity of the metal. We take 
p = nA~^. The pressure can now be expressed in terms 
of the equilibrium pressure Pq, and the perturbation to 
the electron number density per unit length, 6n = SpA^, 
as P = Po + A~'^MB^Sn. Now B is of order of the Fermi 
velocity of the metal composing a wire. The pressure 
force per electron, to be added to the right side of ([3]) , is 
then 

1 dP ^(fR(X,t) 

F = = MB^ (8) 

podX dX^ ^ ' 

where from the linearized continuity (Eq (H])), the pertur- 
bation to the electron number density has again been ex- 
pressed as a function of the equilibrium number density, 
no and the displacement, R of an electron fluid element. 



C. Dispersion Interaction for Insulating wires 

The prior analysis, applying a coupled plasmon ap- 
proach, accounts for conduction of electrons and so pro- 
vides a more complete description in evaluating the vdW 
interaction between two conducting coUinear wires. How- 
ever, we would expect the pairwise summation method 
of individual R~^ atomic contributions, which does not 
allow for electron movement along the wires, to remain 
an apt description of the dispersion interaction between 
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FIG. 2: Pairwise summation of an element of wire A with all 
corresponding elements on wire B 



insulating wires with a sufRciently large gap. We there- 
fore applied two different modifications to the theory de- 
scribed above: each is designed to explore the difference 
between metallic and non-metallic wires. 

(A) A modified pairwise additive approach was ap- 
plied to calculate the dispersion interaction between two 
coUinear insulating wires. Using (O, ^ and 0, we 
first calculated numerically and tabulated the dispersion 
interaction energy Ey(m'{£, d) of short metallic wires of 
length L — £A with separation D = dA between the 
ends. (For fixed £ « d we found Eydw{^,d) to be of 
the conventional form —C{£)d~^ as expected: see the Re- 
sults section below). We then modelled a long insulating 
wire of length L = i£A, {i = 1,2,3,...) by a collection 
of short wire segments of equal length £A, placed end to 
end, but with electrons not allowed to flow from segment 
to segment. This crudely represents individual atoms 
or bonds on an insulating wire, as electron movement 
along the wire is restricted, now being confined to within 
each wire segment. Assuming the mutual polarization 
between any two wire segments is not affected by any 
other wire segment composing the collinear wire system, 
the pre-tabulated dispersion interaction E^dwi^^ d = q£) 
can be used between any two segments of length £A lying 
in opposite wires, and separated by distance q£. Summa- 
tion of the dispersion interaction between each segment 
on one wire and all segments on the opposite wire (as il- 
lustrated in Fig. [2|) then gives the total vdW interaction 
for insulating wires, separated by distance D — p£, in 
the pairwise summation approach: 



-^J^R{X)^ I R{X')x 



"oe" jd/2 

-l\2 I /12 



-2{X-X'Y+A' ^ -2[X + X'f + A^ \ . 

i TTTT dX 



^{X - X'Y + A^f'^ {{X^X'Y+A^) 
MB^ d^R 



5/2 



noe^ dX"^ 



Each eigenfrequency of (fTO|l is of the form 



" \ / " "i, I " pin 



(10) 



(11) 



where Cli is an eigenfrequency of (jlOp with the pinning 
term absent. A semiconducting or insulating wire with 
an electronic energy gap Eg can be modelled by includ- 
ing a finite pinning frequency rJpm 7^ 0. A metallic wire 



is modelled by setting Vlj, 



0. The idea is that each 



electron experiences a force tending to tie it to its la- 
belled equilibrium position X. In this simple approach 
the pinning energy hftpm is loosely associated with the 
electronic gap Eg of an insulator, Mlpin « Eg. The 
pinning has the effect of reducing the polarizability of 
the electron gas, stiffening the plasmon modes and mak- 
ing them less affected by coupling between the wires. 
The modified frequencies ITTI) are then substituted into 
^ to obtain the vdW energy as before. A similar 
pinning approach gives the conventional power laws for 
the vdW attraction between insulating two-dimensional 
sheets [E^dw oc — _D~^) and for distant parallel insulat- 
ing wires^ {Eydw oc —D^^) . As we will show in the 
following Results section, for the present problem in the 
appropriate parameter regime the pinned approach also 
produced results very similar to approach (A) in which 
the wire was cut into short segments. Both will be seen 
to give results for the van der Waals interaction between 
collinear wires that is very different from the case of 
metallic wires. 



EZw""''"iD=p£)= J2 E,dwii,d={p + m + n)£). 

m.n— 

(9) 

Here each wire has length L = i£ where i is an integer. 

This method provides an increasingly accurate repre- 
sentation of two collinear wire insulators with decreasing 
length of the individual wire segments, as electron move- 
ment is restricted to the individual wire pieces. 

(B) We also explored a more conventional approach to 
model an insulator, by adding a harmonic pinning force 
— Afr2pj„i?(X) to the right-hand side of the equation of 

motion for MR. This, plus transposition and inclu- 
sion of the pressure term ([5]) resulted in an eigenvalue 
problem of the form (cf ^) 



III. NUMERICS 

To solve Eq (fTUl) numerically, the integral is discretized 
by application of the trapezoidal rule, which along with 
a zero boundary condition for the displacement at the 
edges of each wire, converts the integral equation (fTO|) 
to a symmetric matrix eigenvalue equation of dimension 

— 1, where N is the total number of partitions of the 
integral term in Eq PH]) . With inclusion of a pinning 
force, the matrix eigenvalue problem is of the form 



M(»2 - n^^ 



N-l 



(12) 
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for 1 < i < — 1. Here Gim is the matrix 

_ ~2{{i-m)L/Nf 
(((^-m)L/7V)2+A2)5/2 

^ -2{D+{i + m)L/Nf 
((i? + (i + m)L/iV)2 + A2)5/2 ■ 

Here = L/N, 1 < m < - 1, is the weight- 

ing function for trapezoidal integration: the end-weights 
Wq — Wn — L/2N are not sampled because of the pin- 
ning condition Rq = Rn — 0, so that the {N—1) x (N—l) 
matrix equation (|12p is symmetric. 
In (Ha) 

Pirn — ^i—l,ra '^^i,m ~^ ^i+l,m (-i*^) 

is a discrete dimensionless version of the second- 
derivative operator. In obtaining (|12p from piT)) . X and 
X' have been replaced by discrete positions Xi and Xm 
defined in terms of the absolute wire length L, wire sep- 
aration distance D and N, by Xi = D/2 + iL/N for i — 
0, 1, 2....N and X„ = D/2 + raL/N for to = 0, 1, 2....N. 

A. Dimensionless form of equations 



is a dimensionless measure of the importance of the pres- 
sure term. 

Equation ^ holds for 1 < i < iV - 1. 

The discrete eigenfrequencies VLj from numerical solu- 
tion of the {N — 1) X [N ~~ 1) matrix equation (15) were 
used as follows in calculating the attractive dispersion 
energy E^dw between the wires (c.f. (7)): 



t JV — 1 

E.dw = 2 5^ + nY\d) ~ 2nf\d -> oo)) . 

(18) 

Here VL^j^\d) is the j*'' eigenvalue of (15) with the -|- 
sign used in (16), corresponding to even plasmon modes. 
Similarly fjj ^ (d) is the eigenvalue with the — sign in 
(15), corresponding to odd modes. fl^^\d oo) is the 
jth eigenvalue with only the first term kept in (16), cor- 
responding to infinite separation d —^ oo. Note that the 
even and odd modes are degenerate in the limit d ^ oo, 
acounting for the factor 2 in the last term of (fT8|) . 



The eigenvalue problem can be written in dimension- 
less form as follows. Firstly we scale all lengths by the 
effective radius A, introducing dimensionless quantities 
d, £, Ti and Wi such that D — dA, L = £A, Ri = ViA 
and Wi — WiA. We also define Jim = A^dm which is 
dimensionless. We introduce a characteristic frequency 
flo defined by 



MA^ 



(14) 



which is \^p, where ^Ip = ^ A'n:nQ{'K A^)~'^ / M is the 
bulk plasma frequency calculated assuming a bulk elec- 
tron density n^o = no{'KA^)~^, as though the ID density 
no were distributed uniformly throughout a cylinder of 
radius A. Then the equation of motion \12^ can be writ- 
ten in dimensionless form as 



pin 



N-1 ^ 
m— 1 ^ 



2^ 



(15) 



where 



Jii 



± 



1 - 2^2(^ _ mfN-"^ 

(£2(i-TO)27V-2-H 1)^/2 

l-2{d + t{i^m)N-'^f 
572 



(16) 



{{d + l{i-'m)N-^f + 1^ 

for even (+) and odd (-) solutions respectively. In (UHl), 

52 



/3^ = 



A^nl 



(17) 



B. Numerical Convergence 



Numerical convergence of the dispersion interaction 
with respect to N was also verified. A larger N corre- 
sponds to dividing the integral into more partitions; con- 
sequently, as more "points" along the wire are sampled, 
larger N means plasmon modes with shorter wavelengths 
are also sampled. The lowest plasmon mode wavelengths 
are expected to contribute negligibly to the dispersion 
interaction, so for sufficiently large N numerical conver- 
gence of the vdW interaction should be evident. We 
found N up to 800 was sufficient for convergence of the 
dispersion interaction for wires separated hy D = 2A at 
all of the lengths treated here (up to 80^). Increasing the 
internal pressure, corresponding to using a larger Fermi 
velocity, gives a more rapid convergence of the dispersion 
interaction with respect to N and a weaker vdW attrac- 
tion between the two wires. This is because the pressure 
term stiffens the short-wavelength modes, making them 
less sensitive to the weak interaction between the wires. 
We also found that inclusion of the pressure improved 
the numerical stability. 

The value of N necessary for convergence (to a given 
fractional level) increases with wire length and also with 
wire separation distance. For comparison of the disper- 
sion energy for conducting and insulating wires, and in 
investigating variation of separation distance for wires 
much longer than the separation distance, we employed 
N sufficiently large to achieve convergence of the vdW 
interaction to about three significant figures. 



7 



Choice of numerical parameters for carbon 
nanotubes. 



The basic approach described above should apply to 
any system where there is an essentially one-dimensional 
electron gas with long mean free path, such as a metallic 
nanowirc or the conduction band of a metallic (n,n) car- 
bon nanotube. The results using (jlSp are dimensionless 
and therefore universal in a sense. However the disper- 
sion energy from ([7]) is thereby obtained in units of h Hq 
defined in (fT4)) . and this quantity will be different for 
different ID systems. Our dimensionless lengths are de- 
fined in terms of the Coulomb smearing length A, which 
is taken to be the radius of the quasi-lD electron gas, 
e.g. the tube radius in the nanotube case, ft is also nec- 
essary to choose the dimensionless input parameters /?, 



Opj„/fio to suit the particular system. 



For the numerical work below we chose parameters 
roughly appropriate to the conduction band of a (5,5) 
single-walled carbon nanotube. While a metallic nan- 
otube is an excellent metal, it is not a free-electron metal. 
We therefore first establish a correspondence between our 
classical Newton II approach and the relevant g — > 
quantal response by considering an infinitely long wire 
as follows. A simple model of the rigid time-oscillatory 
displacement of the k-space Bloch electron distribution 
under an external electric field Eq exp^iqz — icut) gives 
the bare density-density response for g ^ as 



(19) 



Here the Fermi velocity vp = fT''^^de{k)/dk\k=kF is the 
gradient of the 1-dimensional Bloch energy e{k) at the 
Fermi point k ~ kp. vp can be calculated by differentia- 
tion of the analytic Bloch electron dispersion for single- 
walled (n,n) nanotubes given in Eq (4.6) of Saito and 
Dresselhaus^. In this way we obtained the Fermi ve- 
locity of the conduction bands (those with q = n and 
g — 2n in the notation of Saito and Dresselhaus) as 
Vp — \/3a |t| /(2fi.). Here \t\ — 3.03 eV is the hopping 
parameter and a — 4.65 Bohr radii is the length of the 
primitive translation vector of the parent 2D graphene 
lattice. This gives vp = 9.9 x lO^m/s independent of 
n, a value quite comparable to Fermi velocities in 3D 
metals. 

The bare response of the infinite wire can also be calcu- 
lated via the Newton II approach (used in our numerical 
treatment outlined above for the case of finite length), 
giving 



Xo{q,^) = noq^ /{muj'^) 



(20) 



where no is the ID density of the electron gas and m is 
the effective classical mass. 

By comparing (I^Dl) and (fTOl) we establish the equiva- 
lence 

rLQ/m = vp/iirh) . 




FIG. 3: Dispersion energy (f^^E"'''^ vs. d for dimensionless 
length £ — 2, pressure parameter /3 = 1.38, and zero pinning. 



Thus the energy fi,r2o that scales our vdW energy predic- 
tions (see Eqs (HH), ^ and dTS])) is 



I VP e2 8.02 



(21) 



where we have used the result''^ A — \/3na/{27r) for the 
radius of a (n,n) carbon nanotube. 

Assuming that the characteristic velocity B entering 
the pressure term Eq ([8]) is of order vp, we find that the 
dimensionless pressure parameter (3 appearing in ()15p is 
of order unity. We used the value (3 = 1.38 throughout. 

We also need to choose values for the pinning frequency 
mimic the effect of an insulating energy gap. 
Semiconducting carbon nanotubes have electron energy 
gaps up to 0(1 eV). Assuming that Mlpm is of this 
order, we obtain from the numbers above a ratio LOpin = 
^pin/^o = 0{n/8) for a (n,n) nanotube. We explored 
cases with values of Wpin varying from (representing the 
conduction electrons of a metallic wire) to = 1.73 
(representing strongly semiconducting electron bands). 



IV. RESULTS AND DISCUSSION 

We first verified that, in the "pointing" geometry 
of interest here, our numerical coupled-plasmon ap- 
proach yields the expected dispersion interaction of form 
Evdw ~ —Ce{£)d~^ for two wires whose (equal) dimen- 
sionless length i is much less than their dimensionless 
separation d measured between their ends. In this limit 
the d~^ law can be predicted analytically by modelling 
each piece of wire as a polarizable dipole^i^. To test this 
power law explicitly, in Fig. [3] we plot d^^Eydw vs d, 
where c?cm = d -I- ^ is the distance between the centers of 
the wires. 

For each chosen value of £, the plot settles down to a 
constant value —Ce{£) for d >> £, which confirms the ex- 
pectation. This basic behavior is observed both for 
metallic wires and for semiconducting wires (modelled 
by adding a pinning force as described above). For a 
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FIG. 4: Dispersion energy E^dw vs. wire length I for fixed 
separation d = 2 and pressure parameter (5 = 1.38 
Diamonds: metallic wire {ijjpin ~ 0). Solid squares: Semi- 
conducting wire with Upin — 0.707. Solid triangles: semi- 
conducting wire with Upin = 1.0. Crosses: semiconducting 
wire with ujpi„ = 1.224 . Asterisks: semiconducting wire with 
uipin = 1.732. Solid circles: Metallic wire (cUpin = 0) that has 
been cut into pieces of length 2. The unit of energy (on the 
vertical axis) is hQo as defined in Eq (|14|l and (|2ip . For the 
conduction band of a (5,5) carbon nanotube, this energy unit 
is 1.6 eV. 



given wire length £, the coefficient Ce{£) depends consid- 
erably on the pinning frequency Upm i.e. metallic wires 
i'-^pin = ) attract more strongly than insulating ones 
with the same concentration of van der Waals interact- 
ing electrons, but still with a d^^ law for the energy. 

Much more interesting is the behavior of the disper- 
sion energy when the separation d is less than the length 
£. Here a single dipole does not adequately represent the 
electronic response of a wire, and consequently no single 
power law d^ emerges for the energy E^dw{d) because a 
multipolar expansion would be required in order to gen- 
eralize the dipolar argument. We are therefore out of the 
"asymptotic" region of a distant vdW interaction, though 
still not in a region of electronic cloud overlap. In this 
regime, a conventional approach would be to sum (inte- 
grate) contributions of form —C\ri — ^ from small 
segments of the two wires at positions ri and r2 respec- 
tively. This approach predicts that, for fixed d, the dis- 
persion energy for increasing wire length i will saturate 
rapidly as soon as £ exceeds d, to a value proportional to 
—d~^. We find that this saturation does occur for insulat- 
ing wires, as modelled either with a substantial pinning 
frequency Wpi„, or by cutting a metallic wire into mutu- 
ally insulated pieces shorter than the separation. It is 
definitely NOT true for intact metallic wires {tOpm — 0) 
as modelled by our full numerical coupled-plasmon ap- 
proach. In the metallic case the interaction continues to 
grow with wire length and still has not saturated for wires 
many times longer than the separation. These findings 
are illustrated in Fig. |4l 

We also investigated the behavior of our coupled- 
plasmon interaction right down to zero separation of the 
wires: see Fig. [S) Of course a realistic treatment of this 
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d 

FIG. 5: Coupled-plasmon interaction at small wire separation 
d. 



regime would require the inclusion of covalent and other 
forces that are not described in our formalism and that 
are specific to the particular quasi- ID system. Further- 
more a detailed plasmon treatment even of the dispersion 
part of the energy in this regime would require electronic 
response functions beyond the long-wavelength descrip- 
tion that we have used. However our results establish the 
important result that our method saturates naturally to 
a finite value at contact, unlike empirical van der Waals 
correction schemes of the ^ CijR~j^ type, which have to 
be cut off in a somewhat arbitrary manner to avoid di- 
vergence. Another significant conclusion from our num- 
bers in this small-separation regime concerns the relative 
strength of the dispersion energy for metallic and non- 
metallic cases. Because the analytic results in various 
geometriesi^i^ for the metallic interaction fall off with 
a lower power and hence dominate the insulating result 
at large distance, one might speculate that the opposite 
ordering holds at small separation - i.e. that the insulat- 
ing interaction is stronger than the metallic one at small 
separation. However such an argument depends on the 
assumption of the same single power law for small sep- 
arations as well as large. This is not the case, and at 
small separations d a power series rather than a single 
power would be required for this type of analysis. In fact 
our results show that the metallic attraction is stronger 
than the insulating one at all separations, for constant 
electron number density. 



V. SUMMARY AND CONCLUSIONS 

We have explored the dispersion interaction be- 
tween two collinear one-dimensional electronic systems 
("wires") separated by a vacuum gap, in the elec- 
tromagnetically non-retarded regime. We considered 
the metallic case and also two models of an insulat- 
ing/semiconducting case. The metallic case was mod- 
elled via free one-dimensional electrons with an electronic 
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pressure term, but this model was also matched to a 
quantum band model for the electrons for the purpose of 
modelling the conduction electron bands in carbon nan- 
otubes. For the insulating cases, either a harmonic pin- 
ning force was added to the equations, or a model of very 
short, mutually insulated metallic pieces was used, with 
rather similar results to the case of a pinning force. Our 
calculations provide significant insight additional to that 
previously obtained from asymptotic dispersion energy 
calculationsi^iii^i^i^i^ on low-dimensional systems with 
at least one infinite spatial dimension. These previous 
analytic calculations were all performed for the asymp- 
totic regime where a single power law (sometimes multi- 
plied by a logarithm) suffices to describe the dispersion 
energy as a function of separation D. These previous 
works showed that the distant interaction in these highly 
anisotropic systems shows major qualitative differences 
between the metallic and insulating cases. Specifically, 
the power law is different in the two cases, with the metal- 
lic interaction being stronger in the asymptotic regime 
studied. 

In contrast to those fully asymptotic calculations, we 
treated wires of finite length L, with a finite separation 
D between their ends, a regime where a single power 
law in D does not describe the interaction. Unlike the 
asymptotic cases, we find that the interaction at finite 
separations is a continuous function of the pinning fre- 
quency (i.e. a continuous function of the electronic en- 
ergy gap). Our coupled- plasmon dispersion energy does 



however remain a strong function of the gap right down 
to the limit of contacting wires, where it remains finite, 
unlike asymptotic formulae for the same quantity. In par- 
ticular we find a large difference between the attractive 
dispersion energy of the most polarizable electron bands 
of metallic and semiconducting carbon nanotubes. The 
dispersion energy was found to decrease monotonically 
with energy gap at all values of L and D that we stud- 
ied. 

We also found that a conventional pairwise summation 
of terms of form —CR~^ is not an adequate description 
of metallic wires, with a major discrepancy between this 
model and our coupled-plasmon results. Specifically, the 
coupled-plasmon energy, while similar to the summation 
approach for very short wires, continues to grow with 
length L at fixed separation _D, while the — SCi?"^ en- 
ergy saturates rapidly to an L — > oo value that is much 
smaller than the coupled-plasmon energy. 

A more detailed microscopic model, including a de- 
tailed orbital description with Pauli repulsion and non- 
dispersive bonding forces, will be required to describe the 
cohesive forces at distances where the electronic clouds 
overlap. The present results suggest, however, that long- 
wavelength coupled charge fluctuations will still enhance 
the van der Waals force at these distances, necessitating 
a highly nonlocal model of the Coulomb screening even 
in this limit. Theories that yield a ^ R"^ asymptotics 
cannot be relied upon to achieve this. 
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